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The virial theorem for nonlinear problems 2 

Abstract. We show that the virial theorem provides a useful simple tool for 
approximating nonlinear problems. In particular we consider conservative nonlinear 
oscillators and a bifurcation problem. In the former case we obtain the same main 
result derived earlier from the expansion in Chcbyshev polynomials. 



1. Introduction 

In a recent paper Belendez et al [1] showed that the widely used small-amplitude 
approximation cannot always be successfully applied to nonlinear oscillators. To 
overcome this difficulty the authors proposed the expansion of the nonlinear force in 
terms of Chebyshev polynomials. This alternative linearization of nonlinear problems 
proved to be remarkably more accurate and efficient than the straightforward small- 
amplitude approach. Besides, the Chebyshev series applies even to such difficult cases 
where the Taylor series fails [1]. 

The purpose of present article is to discuss an alternative approach to nonlinear 
problems based on the well-known virial theorem [2]. In Sec [2] we outline the main 
results of Belendez et al [1] for conservative nonlinear oscillators. In Sec. [3] we develop 
the virial theorem, apply it to conservative nonlinear oscillators, and compare its results 
with those obtained by Belendez et al [1]. In Sec. H] we apply the virial theorem to a 
nonlinear problem that exhibits bifurcation and compare its results with the exact ones 
and with those produced by the small-amplitude approximation. In Sec. |5] we discuss 
the main results of the paper and draw conclusions. 

2. Conservative nonlinear oscillators 



Belendez et al [1] considered nonlinear conservative autonomous systems given by the 
second-order differential equation 

x + /(x) = (1) 



The virial theorem for nonlinear problems 3 
with the boundary conditions x(0) = A, x(0) = 0. Here a point indicates derivative 
with respect to t. In particular, Belendez et al [1] restricted themselves to odd functions 
/(— x) = —f{x) that satisfy xf > 0. 

The approach proposed by Belendez et al [1] consists in the expansion of the force 
in a series of Chebyshev polynomials of the first kind T„(z): 

oo 

/(x) = EWl(^mn+l(y) (2) 
n=0 

where y = x/A. These polynomials are given by the recurrence relation 

To(^) = 1 
T,{z) =z 

Tn+l{z) = 2zTn{z) - Tn-l{z) (3) 

and are orthogonal in — 1 < 2; < 1 with the weight function w{z) = (1 — z'^)~^^'^: 

fa - Z^)-^'^T^{z)Tr,{z) = |(1 + 5^,)5mn (4) 

Therefore, the coefficients of the expansion ([2]) are given by 

fe2n+i(A) = -/'(!- y')-"'T,^n+Mf{Ay) dy (5) 

Notice that there is a misprint in the weight function shown by Belendez et al [1]. 

If we keep only the first term in the expansion ([2]) the differential equation ([T]) 
becomes that for a harmonic oscillator 

i + ^^ = (6) 

with a frequency 



that depends on the amplitude A. This expression proves to be remarkably accurate for 
many problems [1] in spite of its simplicity. 
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Here we consider the same differential equation ([T]) with the more general boundary 
conditions 

x{h)x{h) - x{a)x{a) = (8) 
If we integrate the equation 

dt •' ^ ^ 



we obtain 



nl x''-^x^dt= j\'^fdt + x{hYx{h)-x{aYx{a) (10) 



In particular, when n = 1 we have 

b fb 

x^dt= / xfdt (11) 

because of the boundary conditions ([8]). 

We now apply this general expression to the oscillators studied by Belendez et al [1] 
that are periodic of period r. In this case the kinetic energy is 

(12) 

and if we choose a = and b = t equation fllip becomes the well-known virial 
theorem [2] 

2K = If (13) 
where the expectation values are defined as 

F = - r Fdt (14) 



r JO 

The virial theorem is known from long ago [2]; its name comes from the fact that xf 
is known as the virial of the forces in the mechanical system. This theorem reveals the 
balance between the kinetic and potential energies [2]. 

The exact trajectory x{t) satisfies equation f|T3l) . If we propose an approximate 
trajectory of the form 

Xappif) = Acosijjjt) (15) 
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where uj = 2tt/t is the frequency of the oscillator, then it is reasonable to set this 
approximate frequency so that Xappit) satisfies the virial theorem (fT3ll . If we substitute 
equation (fTS!) into equation (fT3!) we obtain 

..A^ = 2r'\fit='-^f' y^^iy (16) 

Jo UJ J-1 vi — y 

by means of the change of variables y = cos{ujt). This is exactly the equation for the 
frequency ([7j) derived by Belendez et al [1]. 

We appreciate that both the virial theorem and the first term of the Chebyshev 
expansion lead to the same approximate frequency. 

4. Bifurcation 

Equation (fTTj) is sufficiently general for the treatment of a wide variety of interesting 
nonlinear problems of the form ([T]). In this section we consider the Bratu equation 

u"{x) + Ae"(") = 0, m(0) = = (17) 

that appears in simple models for the study spontaneous explosion due to internal 
heating in combustible materials [3,4]. It is also interesting for another reason: it 
is a simple strongly nonlinear problem that can be exactly solved. Therefore, it is not 
surprising that it has become a useful benchmark for testing approximate methods [4-8]. 
It is well-known that the solution to the Bratu equation is [7] 

where ^ is a root of 

cosh(^^/2)2 ^^^^ 
This equation exhibits two solutions when A < Ac, only one when A = Ac, and none 
when A > Ac, where the critical A-value Ac is the maximum of A(^). We easily obtain it 
from the root of dX{6)/d6 = that is given by 

e^=(0c-2)-0c-2 = O (20) 

The exact critical parameters are 6c = 2.399357280 and Ac = 3.513830719. 
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The slope at origin 

displays a bifurcation diagram as a function of A as shown in Fig. [1] (it is not difficult to 
obtain it by means of a parametric plot using equations f|T9|) and fl2Tl) ). For the critical 
value of A we have u'(0)c = 4. 

In what follows we show that the virial theorem is suitable for estimating the 
form of this bifurcation diagram. We simply have to introduce a trial function 
which satisfies the appropriate boundary conditions, into the expression for the "virial 
theorem" 

dx + \ C we" dx = Q (22) 
Jo Jo 

Notice that the exact solution satisfies u"{x) < for all < x < 1; therefore u{x) is 
positive and do not have zeros between the end points. This conclusion will guide us 
towards the choice of the trial function. 

One of the simplest functions that meets the criteria just indicated is 

u{x) = Ax{l - x) (23) 

A straightforward calculation shows that 

4^5/2 

^ = ^^^^ ^ (24) 

3 [V^iA - 2)e^/4erf (v^/2) + 2Va\ 

and the slope at origin is u'{0) = A, so that we can easily plot m'(0) vs A parametrically. 
Fig. [1] shows that this expression is suitable fo the lower branch (small A) but it 
is not so accurate for the upper one (large A). However, it provides a reasonable 
description of the bifurcation diagram and the critical parameters Ac = 3.569086042 
and u'{0)c = 4.727715383 are remarkably close to the exact ones. 

Another simple variational function that meets the required criteria is 

u{x) = Asin(7rx) (25) 

that leads to 

Att^ 

^ ^ 2{2 + 7r[/i(A) + Li(A)]} ^^^^ 
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where Iu{,z) and Li,{z) stand for the modified Bessel and Struve functions [9], 
respectively. In this case ^^'(0) = tiA and Fig. [T] shows that this expression is 
shghtly less accurate than the preceding one for the lower branch and certainly more 
accurate for the upper one. Besides, this trial function yields better critical parameters: 



The Bratu equation is also suitable for revealing the limitation of the linearization 
by means of an expansion in a Taylor series. If we neglect the nonlinear terms in the 
expansion: e" = 1 + m + . . . then we can solve the resulting differential equation exactly 
and obtain 



Fig. [T] shows that this approach based on the Taylor expansion is unable to reproduce 
the upper branch of the bifurcation diagram. The explanation is quite simple: the 
solution for the lower branch is considerably smaller than the one for the upper branch. 
Therefore, an expansion based on small values of u will necessarily produce the former 
and fail on the latter. On the other hand, an expansion in appropriate orthogonal 
polynomials (or the virial theorem) provides an acceptable description of both branches 
of the bifurcation diagram. 

5. Conclusions 

We have shown that the approach derived by Belendez et al [1] from the first term of 
the expansion in Chebyshev polynomials can also be obtained by means of the virial 
theorem. It is clear that we can introduce the approximation in two different ways: 
as the first term of a systematic numerical method or as the requirement posed by 
the virial theorem with a direct physical interpretation. One or the other point of 
view (or perhaps one after the other) may be useful for teaching an undergraduate 



Ae = 3.509329130 and u'(0)c = 3.756549365. 




(27) 



In this case the slope at origin is 




(28) 
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course on classical mechanics. One can easily derive and discuss the virial theorem 
for mechanical problems and then generalize it for the treatment of arbitrary ordinary 
nonlinear differential equations. One advantage of the approach based on the virial 
theorem is that it is also suitable for the treatment of quantum-mechanical problems 
as well [10]. 

The virial theorem provides us with a quite general expression that may be useful in 
the study of many nonlinear problems. As an example we have shown that the approach 
is suitable for the treatment of the well-known Bratu equation that appears in simple 
models for heat combustion [3-8]. In this case we have been able to try two different 
approximate solutions which may probably be more difficult if one merely resorts to an 

expansion in orthogonal polynomials. 
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Figure 1. Bifurcation diagram for the slope at origin u'{Q) in terms of A obtained 
by means of the exact expression (sohd hne), ^(a;) = Ax{l — x) (dashed hne), 
u{x) = ^sin(7ra;) (dots) and Taylor linearization (circles) 



